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The vortex core radius p„, defined as the peak position of the supercurrent around the vortex, 
has been determined by muon spin rotation measurements in the mixed state of La2-a;Sra;Cu04 for 
X — 0.13, 0.15, and 0.19. At lower doping {x = 0.13 and 0.15), pv{T) increases with decreasing 
temperature T, which is opposite to the behavior predicted by the conventional theory. Moreover, 
pv{T — > 0) is significantly larger than the Ginsburg-Landau coherence length determined by the 
upper critical field, and shows a clear tendency to decrease with increasing the doping x. These fea- 
tures can be qualitatively reproduced in a microscopic model involving antiferromagnetic electronic 
correlations. 



The structural and electronic properties of the flux line 
lattice (FLL) state in high- Tc cuprates have been exten- 
sively studied by microscopic techniques, including scan- 
ning tunneUng spectroscopy (STS), small angle neutron 
scattering (SANS), nuclear magnetic resonance (NMR), 
and muon spin rotation (/iSR) . Rich information has been 
provided not only on the superconductivity itself, but 
also on the unique FLL state as a 'vortex matter' real- 
ized in YBa2Cu307_5 (YBCO) and BizSraCaCuaOg+s 
(BSCCO) systems. In particular, recent developments 
concerning the /iSR experimental technique have made it 
feasible to obtain microscopic details of the spatial field 
distribution B{y) in the FLL state directly from the /xSR 
time spectral. The application of this advanced /xSR 
technique to various type II superconductors including 
YBCO has revealed that the vortex core radius py, which 
must be proportional to the Ginzburg-Landau (GL) co- 
herence length near the transition temperature {Tc), 
strongly depends on the external magnetic field: p„ ex- 
hibits a sharp increase with decreasing field at lower 
fields, while it tends to converge to ^gl = v^<f>o/27riJc2 
at higher fields (with $0 and Hc2 being the fiux quan- 
tum and the upper critical field, respectively)^^. This 
indicates that the conventional picture of the FLL state, 
in which the vortices are regarded as being simple arrays 
of rigid cylinders (with a radius ^gl) containing normal 
electrons, is valid only near Tc, i.e., the region prerequi- 
site for the application of the GL theory. A recent calcu- 
lation based on the quasiclassical Eilenberger theory^ has 
successfully reproduced the observed field dependence of 
Py in GeRu2 in which an isotropic order parameter is re- 
alized for H/Hc2 < 0.5^. On the other hand, the result of 
a similar calculation for anisotropic d-wave superconduc- 
tors is far from satisfactory to explain the experimental 
observation in YBCO, where py exhibits a much steeper 
change compared with the theoretical prediction. 

Meanwhile, similar studies on the La2-a;Sra;Cu04 



(LSCO) system are still in a preliminary stage because 
both the STS and SANS techniques have rarely been suc- 
cessful in observing vortices in this compound. It is only 
very recently that the SANS signal from FLL has been 
reported in an overdoped LSCO with the revelation of 
a square FLL structure at higher magnetic fields^. The 
previous /xSR studies on LSCO were mostly concerned 
with remnant magnetism over the underdoped region, in- 
cluding that near x = 1/8^, or with magnetic penetration 
depth. A, which was evaluated from the total lincwidth 
without any detailed modeling of the microscopic struc- 
ture of vortices^. To our knowledge, there has only been 
one attempt to derive the vortex core size, which found a 
relatively large cutoff parameter assuming the modified 
London model^ . 

In this work, we report on high-precision /iSR mea- 
surements of the magnetic field distribution in the mixed 
state of La2-a;Sra;Cu04 single crystals. Wc found that 
the field distribution at lower fields {H < 0.4 T) is per- 
fectly reproduced by a modified London model with a 
Lorentzian cutoff. The actual size of py defined as the 
peak of the supercurrent density deduced from Maxwell's 
relation was quite large (e.g., pyC^ 80 A for x = 0.15), be- 
ing much larger than ^gl (— 23 A for x — 0.15). More 
importantly, we have found that i) py{T) in LSCO with 
X = 0.13 and 0.15 tends to increase with decreasing tem- 
perature T, and that ii) pv{0) exhibits a monotonic de- 
crease with increasing doping x. While these features are 
qualitatively unique compared with YBCO^" and have 
no simple explanation within the conventional theoreti- 
cal model, they can be reproduced in a model involving 
strong antiferromagnetic correlations competing with su- 
perconductivity. 

/iSR experiments on two sets of single crystalline 
La2_a;Sra;Cu04 (one with x = 0.15, labeled LA15, and 
those with x =0.13, 0.15, and 0.19, labeled LB13, LB15, 
and LB19, respectively) were performed on the M15 
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FIG. 1: a) Time differential fi-e decay asymmetry in 
Lai.85Sro.i5Cu04 (LA15) at 15 K and = 0.2 T, where 
the solid curve shows the result of a fitting analysis with a 
Lorentzian cutoff (see text), b) Difference between the data 
and the solid curve, which is defined as the error of the fitting 
analysis. 

beamline at the TRIUMF muon facility which provides 
a beam of nearly 100% spin-polarized positive muons of 
momentum 29 MeV/c. Some results on LB15 were pre- 
viously reported^. The specimen was loaded onto a He- 
flow cryostat and cooled from a temperature above Tc 
after setting the magnetic field at every field point (i.e., 
field-cooling) to minimize the effect of flux pinning. The 
initial muon spin polarization was perpendicular to the 
magnetic fleld H, and thus to the FLL in the supercon- 
ducting state, as the crystal c-axis was alligned parallel 
to H so that the fleld distribution associated with the 
FLL was determined by the in-plane penetration depth 
Xab and the coherence length ^ab- In zero field, //SR re- 
sults from the specimen LA15 exhibited slightly enhanced 
spin relaxation below ~10 K, which was virtually absent 
in LB15. The spin relaxation rate in LA15 at 3.5 K was 
greater than that at 41 K by 0.026(2) /is~^, indicating a 
weak remnant magnetism. In addition, we observed an 
additional component undergoing a fast spin relaxation 
in LA15 in measurements with an applied field above 0.4 
T, which became more prominent at higher fields. In 
the present study, we analyzed data obtained only for 
H = 0.2 T which do not exhibit such additional relax- 
ation, in order to avoid any complication due to mag- 
netism. We will report on high-field /iSR measurements 
in all of these specimens in a separate paper^^, and dis- 
cuss the possible relation to field-induced magnetism sug- 
gested by neutron diffraction^^-^'^'^^. 

Since muons stop randomly along the length scale of 
FLL, the time evolution of complex muon polarization 
P(t) provides a random sampling of the internal field 
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FIG. 2: a) Example of the reconstructed field distribution 
{Bz{r), solid line with data shown by open squares) using 
the result of Fig. Q which yields the spectral density {n{B)) 
in b). The radial distribution of the supercurrent density, 
J(r) = |rotB(r)|, is shown in c) along with a definition of 
the core radius pv For a comparison, those corresponding to 
p„=^GL are also displayed. 

distribution, B(r) = (0,0, S(r)), 

P{t) EE P,{t)+iPy{t) 

/OG 
n{B) exp{ijf,Bt - itp)dBll) 
OO 

n{B) = {5{B - B{r)))r, (2) 

where Px^y{t) is proportional to the time-dependent /i"*"- 
decay asymmetry Ax^y{t) deduced from the corre- 
sponding sets of positron counters N^(t) {ip = wk /2, with 
n = 0, 2 and 1, 3) after subtracting constant backgrounds 
estimated from data at earlier times (i < 0), 

Ax{t) = No{t)lN^{t) = A^P,{t), (3) 
Ay{t) = aN^,2{t)/N^^/2{t)^AoPy{t), (4) 

with a = [iVo(0)iV3^/2(0)]/K(0)iV,/2(0)], is the nor- 
malized instrumental asymmetry, dp is the additional re- 
laxation due to random flux pinning, n{B) is the spectral 
density for the internal field defined as a spatial average 
(Or) of the delta function, 7^ is the muon gyromagnetic 
ratio (=27rx 135.53 MHz/T), and V' is the initial phasei^. 
These equations indicate that the real amplitude of the 
Fourier-transformed muon precession signal corresponds 
to n{B) with an appropriate correction of Cp. In the mod- 
ified London model, B{r) is approximated as the sum of 
the magnetic induction from isolated vortices to yield 

B{v) = B,Y. ^^ll F[K,^,), (5) 
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where K are the vortex reciprocal lattice vectors, Bq 
(~ H) is the average internal field, A is the London pene- 
tration depth, and F{K, is a nonlocal correction term 
with being the cutoff parameter for the magnetic field 
distribution; it must be stressed that is a parame- 
ter to describe the electromagnetic response of a vortex, 
thereby not necessarily equivalent to the core radius p^. 
Considering small anisotropy predicted by the theory'*, 
we assumed an isotropic vortex core and associated field 
distribution near the vortex center. 

TABLE L A comparison of the central field between points 
above and below the superconducting transition temperature 
Tc, where the external field was set to 6 T above Tc and 
then field-cooled. The relative precision of the field (|A///|) 
was better than fO"'' for all cases. Note that the effect of 
field inhomogeneity and/or the Knight shift is larger at higher 
fields. Therefore, the above number gives an upper bound for 
the uncertainty of the central field. 



Samples 


/ = 7^Bo/27r (MHz) 


IA///I 




T>Tc 


r < 4 K (< T,) 




LA15 


812.9905(22) 


812.9902(66) 


0.0(9) X 10-5 


LB13 


813.0765(15) 


813.0533(41) 


2.9(5) X 10^5 


LBf5 


812.8405(17) 


812.7686(30) 


8.9(3) X 10-5 


LBf9 


813.0679(12) 


813.0384(19) 


3.6(3) X 10^5 



As summarized in Table the relative change of Bq 
due to field inhomogeneity and/or the Knight shift was 
less than 10"* throughout measurements below 50 K, 
and therefore Bq was fixed to the value determined in 
the normal state (T > Tc) for the analysis of data be- 
low Tc- This well-defined Bq, together with the strongly 
asymmetric feature of n{B) against Bq, allowed us to de- 
duce the physical parameters including A and with- 
out much ambiguity. Moreover, we have found that 
the deduced values of £,y (and thereby pv) were virtu- 
ally independent of the apex angle (6*) of FLL, while 
A showed a considerable dependence on 9. This fea- 
ture can be readily understood by considering the fact 
that the change in the apex angle has least effect on the 
field distribution of single vortex near the center (cor- 
responding to the high-field end of n{B)), while it does 
modify the overlap of field distribution between the vor- 
tices. Based on this robustness of the analysis result 
against 6, we assumed triangular FLL (0 = tt/S) through- 
out the following data analysis for simplicity. We have 
also found that the /zSR spectra in LSCO at lower fields 
are much better reproduced by the Lorentzian cutoff, 
F{K,£^y) = exp{—^/2K£_y), compared with the conven- 
tional Gaussian cutoff, F{K,^y) — exp{~K'^£,y/2). Fig- 
ure n shows a typical example of the measured decay 
asymmetry, Ax{t) (oc Px{t)), observed at 0.2 T in spec- 
imen LA15 together with the result of a fitting analy- 
sis. The reduced chi-square 'X^ /Nf (where Nf denotes 
the number of degrees of freedom) is close to the ideal 
value of ~ 1 when the Lorentzian cutoff is applied to 
Eq. (O, while it becomes much worse as X^/^f — 1-61 



with the Gaussian cutoff. A similar tendency was com- 
monly observed in all other samples. This observation is 
consistent with the result of a theoretical analysis which 
showed that the Lorentzian cutoff is indeed a better ap- 
proximation at the low-field limitiS*ii. Accordingly, all 
of the /iSR time spectra were analyzed by comparing the 
data with the time evolution calculated by Eqs. Q^^© 
using the Lorentzian cutoff, where Aq, ?/), cTp, ^y, and A 
were free parameters while Bq and 9 were always fixed. 

While the fitting analysis of the /iSR data is performed 
entirely in time domain, one can reconstruct B{r) by 
Eq.(|SJl using the physical parameters deduced from the 
fitting analysis. An example of the reconstructed field 
distribution, B{r) = Bz{r), along the radial direction 
from the core center to a saddle point is shown in Fig. 
together with the corresponding n{B) (Fig. and su- 
percurrent density, J(r) = |rotB(r)| (Fig. |3;). As shown 
in Fig. 13;, the core radius defined by J{pv) = Jmax (where 
"/max denotes the maximum of J{r)) is considerably larger 
than the magnetic cutoff parameter indicating the 
need for a special precaution in interpreting ^„ directly 
as the core radius. However, the result of a data anal- 
ysis for various fields/temperatures indicates that py is 
always proportional to ^t,. It has been estimated that 
^GL — 23 A for an optimally doped sample {Hc2 — 62 
T for X = 0.15i^). Provided that py is independent of 
the field and determined by ^GL(i-e., py— 0.6 ^ O.S^cLas 
predicted by theorjii), then ^„ must be smaller than 13 A 
to reproduce the corresponding core radius (see Fig. [21 . 
However, an attempt to fit the data assuming ^^ = 13 A 
completely fails to reproduce the data; the reduced chi- 
square /Nf was 2.76 at its best when other parameters 
(i.e., Aq, ij}, (Tp, and A) were set free to minimize /Nf. 
This is primarily because, as evident in Fig.|23, the peak 
of n{B) shifts significantly to a lower field, which cannot 
be compensated by any other parameters within the cur- 
rent model. Note that Bq is determined with a relative 
precision better than x 10"'* so that the shift of the peak 
in Fig. \^ (~ 0.2 mT) is readily discernible. Thus, it is 
inferred from these results that py is about three times as 
large as ^gl near the lower critical field. This tendency 
is qualitatively in line with the theoretical prediction, al- 
though the magnitude is yet to be explained^. 

The deduced vortex core radius for samples LB13, 
LB15, and LB19, as a function of the temperature, are 
shown in Figs. Oji-c together with the results for LA15. 
It exhibits a slight increase with decreasing temperature 
for X — 0.13 and 0.15, while an opposite tendency is ob- 
served for X = 0.19. No such anomaly has been reported 
for YBCO*°. Note that the conventional theorjiiSiSS pre- 
dicts a behavior opposite to what is actually observed in 
the former cases (see Fig. EJi). The values deduced from 
the data in LA15 are in good agreement with those in 
LB15, indicating that the influence of the remnant mag- 
netism found in LA15 is not significant a,t H = 0.2 T 
and T > 15 K. The penetration depth extrapolated to 
K was 2559(50) A in LA15, 2446(6) A, 2460(4) A, and 
2051(7) A in LB13, LB15, and LB19, respectively The 
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relatively large error for LA15 is due to the influence of 
the remnant magnetism below 10 K. We also note that 
Op was typically 0.4'--^ 0.5 /^s^^ at the lowest temperature 
in all specimens, showing a common tendency of gradual 
decrease with increasing temperature. This is understood 
as the thermal depinning of FLL from random pinning 
centers. 




LSCO 
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FIG. 3: Temperature dependence of the vortex core radius 
at H = 0.2 T for specimens a) LB13, b) LA15 (triangles) 
and LB15, and c) LB19. d) Calculated pv vs T/Tc using 
the model in Ref.^^, where U /t represents the strength of the 
electronic correlation (Tc is the superconducting transition 
temperature) . The solid curves in a)-c) show fitting results by 
an analitical approximation of d) to deduce the extrapolated 
values pv{Q)- 

For understanding the T-dependence of p^, we per- 
formed a model calculation based on Bogoliubov-de 
Gennes theory, where both the d-wave superconductiv- 
ity and the spatially modulated antiferromagnetic (AF) 
spin correlations are simultaneously considered by incor- 
porating a pairing interaction to the standard Hubbard 
model for two-dimensional square lattice^!. As can be 
seen in Fig. [Sji, the core radius at which the supercon- 
ducting order parameter reaches 60 % of the bulk value, 
is strongly enhanced when the electronic correlation is 
present (i.e., U/t > 0, with U and t being the respec- 
tive on-site Coulomb energy and transfer matrix element 
to the nearest neighboring sites in the Hubbard model), 
while it obeys the prediction by Kramer and Pesch when 
U/t = Oi2iSa. The T-dependence of p„ for U/t > quali- 
tatively agrees with those observed for x — 0.13 and 0.15. 
In the model, this feature is due to AF spin correlations 
induced at the core sites, leading to a reduction of the lo- 
cal pairing amplitude and an associated increase of at 
lower temperatures. The absence of quasistatic antiferro- 
magnetism may be explained by the strongly dynamical 
spin fluctuations in the relevant field/temperature range. 

In order to obtain the extrapolated value Pi,(0), the 
data in Figs.l^L-Sc were analyzed by the ^^-minimization 



TABLE II: Parameter values for Eq. © deduced by fitting 
data in Fig. 



Samples 


pv{0) (A) 


Cl 


C2 


C3 


X'/Nf 


LB13 
LB15 
LB19 


120.9 (4.6) 
95.0(7.3) 
73.3(11.5) 


-6.36(3) 
-6.42(5) 
-5.51(6) 


-2.69(4) 
-2.68(6) 
-2.14(9) 


2.65(2) 
2.72(4) 
2.91(4) 


7.5/5 
9.8/7 
18.5/3 



method using an analytical approximation of Fig. 

p„(T) = p„(0){l+ci(T)+C2(exp(-T)-l)+C3(exp(T)-l)} 

(6) 

under a condition ci < C2 < < C3 with T — T/Tc. 
The result of fitting analysis is summarized in Table 
The values of p^ (0) are plotted against the Sr concentra- 
tion X in Fig. 0] (The data of LA15 were not used for 
this analysis because they were available only for T > 15 
K, leading to large uncertainty in the fitting.) As the 
hole doping progresses with increasing x, the core size 
decreases monotonically. We note that almost equivalent 
result was obtained by an independent analysis using a 
linear function p„(T) = /5„(0)(1 + cT) for T < O.bTc. 
We shall compare this behavior with those of other rel- 
evant parameters. Firstly, p«(0) and have no simple 
correlations, since Tc shows a maximum at x = 0.15. Sec- 
ondly, the superconducting gap. A, does not scale with 
p„(0), if A scales with Tc- In contrast, p„(0) in YBCO 
is larger for the lower Tc, thus following the simple scal- 
ing law^°. It has been noticed in LSCO that the gap (or 
pseudo-gap) increases with decreasing hole concentration 
in the optimum-to-underdoped region. In this sense, the 
(pseudo-) gap value and Pt,(0) vary along with each other. 
This behavior is opposite to that expected in a simple 
BCS superconductor, where the coherence length is in- 
versely proportional to the gap value as ^0 = Jivf/t^^, 
where vf denotes the Fermi velocity. In La2-a;Sra;Cu04, 
the Fermi velocity has been known to have little depen- 
dence on a;. It is interesting to note that the model cal- 
culation shown in Fig.OJi appears to simulate the doping 
dependence if we assume that the incipient AF correla- 
tions become weaker {U/t ^ 0) as the doping proceeds. 
A similar result has been reported for the correlation 
length, £^zn{x), in Zn-doped La2-a;Srj;Cu04 over which 
the pairing is suppresseci^. 

Finally, we point out the potential link between our 
results and those of recent neutron scattering in LSCO 
where a field-induced quasistatic antiferromagnetism has 
been suggestecii2iiii4. Although it has not been clearly 
identified as being due to the vortex cores, they found 
that a long-range AF correlation is recovered in the 
mixed state under moderate magnetic fields of a few 
Tesla. A similar situation is also suggested in YBCO^iSi, 
BSCCOS^, and YBa2Cu40a^. The enhanced vortex core 
radius at lower fields in LSCO may be interpreted as a 
precursor of such a quasistatic correlation. Besides our 
model calculation, theories including those based on the 
t-J model2L2& or SO (5) symmetry'^^ also predict the de- 
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FIG. 4: Doping dependence of pv extrapolated to K at = 
0.2 T and that of the superconducting transition temperature 
(Tc) determined by the Meissner effect. The hnes are guides 
for the eye. 



velopment of AF correlations in place of a superconduct- 
ing order parameter when the latter is suppressed by a 



magnetic field or impurity atoms. When the supercon- 
ducting order parameter is suppressed at the center of 
vortices, one would expect the emergence of AF spin cor- 
relations around the vortex cores, which may be strongly 
dynamical at lower fields. 

In summary, we have found in LSCO that the vortex 
core radius at 0.2 T is about three times as large as that 
estimated from Hc2- The core size increases monotoni- 
cally with decreasing Sr doping with an unusual temper- 
ature dependence for x = 0.13 and 0.15. Compared with 
the cases of other cuprates, these features more strongly 
suggest a possibility that p„ is influenced by the two- 
dimensional AF correlations in this system. 
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